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1. Introduction 

1.1 Motivation 

Twisted mass lattice QCD (tmLQCD) offers an efficient mechanism for eliminating unphysical 
zero modes thereby admitting calculations at lighter quark masses [[T]]. At maximal twist one has 
0{a) improvement [§]. Lighter quark mass calculations require improved statistics for hadron mass 
resolution. This fact and the need to disentangle physical states with the same quantum numbers 
may in part be accomplished through the creation of more operators that represent the channel in 
question. The creation of more elaborate operators also allows for the study of hybrid and exotic 
mesons. 

The violation of parity by the twisted mass action causes many operators to be twist angle 
dependent resulting in correlator contamination by opposite parity states when not at maximal 
twist. This suggests seeking a class of twist independent operators from which the physical parity 
of an operator is readily deduced. 

Operators with displaced quarks have been constructed using group theoretical techniques but 
their usage necessarily requires the calculation of quark propagators from multiple sources. (See, 
for example, Refs. [Q, ^ ^.) Consideration of the operators available through gluonic extension 
alone therefore is calculationally expedient. 

1.2 Lattice Symmetry Group 

While parity (P) and charge conjugation (C) may be conserved by lattice actions, the contin- 
uous rotational symmetry of nature is broken and one requires operators adapted to the symmetry 
group of the lattice. For mesons this is the octahedral group O with 24 elements. The direct product 
of the parity, charge conjugation, and octahedral groups is denoted O^''. 

For an operator adapted to the representation A^*-, where {A\,A2,E,T\,T2} is an irre- 
ducible representation (irrep) of O and P,C G {+,—}, one identifies the possible physical states 
J^'^ to which it corresponds using Table [I] which shows the number of copies «^ of irrep A to 
which the continuum SU (2) irrep J subduces[|6|]. An operator transforming as the irrep E^^ could 
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Table 1: Subduction of continuum irreps to those of the lattice symmetry group. 

have spin content (7^'- = 2^ ,4^ ,5^ , . . .). In theory one needs many operators of each A^*- to 
resolve the many physical spins Z^'- in the tower of states to which A^'^ corresponds. 

1.3 Operator Building Blocks 

One may construct zero-momentum operators transforming as irreps of O^^ from the space of 
operators spanned by 

Mj,k,aAt) = Yjfa{^)Uj{x)Uk{x + hU-j{x + j + k)U^u{x + k)Wh{x) , (1.1) 

X 
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where j^k = ±1, ±2, ±3; 7 7^ k and a and b are spinor indices for a total of 24 x 16 = 384 opera- 
tors. See Figure |l| for a diagrammatic representation. Superpositions of such operators for meson 
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Figure 1: Diagrammatic representation of the building block operator, 
spectra analysis have been suggested and used previously in special cases 

[0,ill|l- Evaluation of 

correlators constructed from operators of this form requires only calculation of propagators from a 
single source. Construction of operators transforming as irrep A^*- is facilitated by observing that 
one can effectively consider the transformation of the spinor and link paths independently and then 
combine them via octahedral Clebsch-Gordan coefficients, as shown in Sections 0, ||, and ^. 



2. Spinor Contribution to Group Structure 

The contribution to the group structure due to spinor indices is determined by the 16 bilinears 
xj/Fxj/, where F represents one of sixteen 4x4 matrices, {/, 75, 74, 7475, "y;-, "K/s, CJ4,-, £,-yi^ayjt}. The 
first four bilinears ai^e scalai^s while the last four three-index objects are vectors under the rotation 
group in the Euclidean continuum. 

Parity (P) and charge conjugation (C) of the bilinears are identified via 

^i^^t = (xj/C'^Y ^VA^f = 74 VA 
^^"f = - {Cliff =^v7^f = TO ' 

where C is the matrix implementing charge conjugation. 

To classify the bilinears into iiTcps of O one uses its character table found in Table |^[|6|]. Since 
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Table 2: Character table giving j ((^ ) for octahedral group O. 

the bilinears are scalars and vectors in the continuum it follows that their respective spans are also 
invariant under O. By inspection it is straightforward to show that the character table for the scalar 
and vector bilinears is identical with that of Ai and Ti respectively. The multiplicity w;^ of irrep A 
in the representation J is given by 

ni = -Y,Pa^''K^rX^'\^), (2.2) 
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where is the number of elements in class t, and go is the order of the octahedral group [Q]. 
It follows trivially by this formula that the scalar and vector bilinears form the basis of A\ and T\ 
irreps respectively. Furthermore it may be verified that the vector bilinear component transforms 
as the row for our choice of matrix representation P'^' '(/?). 

Hence the reduction of the spinor structure of our operators is given in Table |3[ where now 
each bilinear may be classified uniquely by its irrep A^'', row A, and irrep multiplicity index a. 
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Table 3: Octahedral symmetries of the spinor bilinears. Red and green entries are twist angle invariant for 
charged and neutral mesons respectively as shown in Section 5.1. 



3. Link Contribution to Group Structure 

Define the gauge link part of our operators in equation ( |l.l| ) by 

Uj,k{x) = Uj{x)Uk{x + j)U-j{x + j + k)U-k{x + k) 



(3.1) 



where j,k = ibl,ib2,ib3; j ^ k. The parity and charge conjugation contributions due to this link 
structure can then be taken into account by defining the PC-adapted superpositions, 
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Uj,k + PU-j.-k + CUkj + PCU k, j 



(3.2) 



where k is defined via k = ix j and we have suppressed an overall normalization factor of 1/2. 
Diagrammatically, Ufj' is shown in Figure ^ For fixed PC the space spanned by Ufj' is invariant 




PC 




k 



Figure 2: PC-symmetrized basis elements for gauge field links. 

under O and will generate a representation U''^ of the group. A basis for the six-dimensional space 
may be found by restricting to {(1,2), (1,3), (2,3), (2, 1), (3, 1), (3,2)} and this will in turn 



give the particular matrix representation r|^^.j^^^ ^-^ (/?). The character table for U^'^ giving 
is shown in Table Q Use of Formula (22) reduces U'^'^ to iiTcps A^'^ of C^'' as shown in Table ^. 
The symmetrized link fields corresponding to this reduction are 



4 



Toward Group Theory Operators for tniLQCD Hadrons 



Robert G. Petry 



E 3C| 8C3 6C4 6C2 
6 2P {P+l)C 



Table 4: Character table for representation U induced by Uf 
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Table 5: Reduction of U^^ to irreps A'^^ of O' 



ur (-) = I ^'-^^^""^ . ^W , (3.3) 



where (m,?i) and /i are chosen so the denominator does not vanish and the accessible A^*' are 
chosen from Table ^ (See Ref. |1C, page 74].) Here r^^'(/?) is our matrix representation for the 

irrep A. The U^'"'^ (x) transform as the row X of iiTcp A^'^ and ai^e uniquely identified by them. 

4. Total Representation Reduction and Operator Construction 

4.1 Reduction of Direct Product Representations 

Having classified the spinor and link components of our operators into irreps of O^^, it re- 
mains to combine them and reduce all possible direct product representations. Parity and charge 
conjugation for the product representations are given by 

P = PfPu C = CfC,, (4.1) 

where / corresponds to the fermionic part and u corresponds to the gauge part. Noting that the 
character of Ai (g) A2 satisfies ^(^^^-^2) (^^'^ = x''^'K^)X^^^\^) and using Formula (|0), one may 



reduce each direct product into irreps of O. As such the 384-dimensional representation generated 
by the space spanned by Mj^k,a,b{i) may be reduced into irreps of O^'^ as shown in Table ^. Note 
that all possible lattice irreps ai^e accessible, not just the A^^ and Tf^ of simpler local operators. 



4.2 Operator Construction Using Clebscli-Gordan Coefficients 

The construction of the lattice-symmetrized operators themselves is accomplished by combin 

' and U£ 



ing the spinor and link operator components F, ' and Uir using the same formulae for parity 



and charge conjugation given in Equation (4.1) and by combining the irreps of O using Clebsch- 



Gordan (C-G) coefficients for Ai (g) A2. The simplicity of the possible direct product representations 
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Table 6: Reduction of representation generated by span of ) into O^^ irreps. 

allows one to use 
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to determine the C-G coefficients. (See Ref. [|^, page 74].) Here jJ., jJ-i, and /I2 are chosen so the 
denominator does not vanish. The lattice-symmetrized operators, (t), are then given by 
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where the allowed iiTcps for the spinor bilinear and link components are determined by Tables || 
and |[ These fix P and C for the operator while the irreps A of O ai^e those allowed by the C-G 
series in Table ^. Dirac indices on F and color indices on U and both indices on the spinors have 
been suppressed. 

5. Twisted Mass Invariant Operators 

5.1 Operators with Definite Parity 

Consider tmLQCD for a quark doublet. A change of basis (by twist angle ft)) gives 

= exp(/ftJ75T3) ^ ' (^-1) 

which leads to mixing among some, but not all, quark bilinears. For example, 

ujid' = uyidco&(0 + uyiysdiin(o (5.2) 

couples to both vector and axial mesons unless CO is carefully tuned, but 

u o^id' = uo^id (5.3) 

couples to the vector and never to the axial. Explicit calculation shows that tmLQCD splits the 
results of Section ^ into two sets according to flavor. Structures in red in Table ^ are only inde- 
pendent of twist angle for chai^ged mesons, while green structures are only independent of twist 
angle for neutral mesons. Despite these restrictions, one finds that all AF^ combinations can still 
be obtained by factoring in some nontrivial link structure. Combining the twist-invariant spinor 
structures from Table || with the link structure in Table ^ results in operators of the irreps outlined 
in Table ^ being accessible. 
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Table 7: Twist invariant irreps for charged and neutral mesons respectively. *lt is to be noted here that 
charged mesons are really only eigenstates of g-parity with eigenvalue G = C( — 1 but we use C for simple 
correspondence with the previous arguments. 



5.2 Preserving Parity in the Correlators 

Building states of definite A^*- might seem insufficient for tmLQCD, since the theory does 
not preserve P. However, tmLQCD does conserve a "twisted parity", P, defined by the product of 
P and {(O —5- —ft)). Therefore, states of definite twisted parity are automatically preserved in any 
correlation function. Since our quark bilinears are independent of (O, P and P are equivalent for 
these bilinears. 
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